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1 .  INTRODUCTION 


In  recent  times  there  has  been  much  interest  in  the  generation  of 
spatial  grids  by  numerical  methods  either  by  algebraic  means  or  by 
solving  certain  partial  differential  equations.  A  collection  of 
various  methods  has  appeared  in  [1]  and  reference  may  also  be  made  to 
the  review  articles  [2],  C3]»  and  to  various  conference  proceedings, 
[4,53  • 

This  paper  is  directed  to  the  problem  of  grid  generation  in  a  given 
surface  by  utilizing  a  set  of  elliptic  equations.  The  proposed  mathe¬ 
matical  model  has  deliberately  been  made  to  depend  on  the  formulae  of 
Gauss,  which  naturally  involve  the  partial  derivatives  of  the  Cartesian 
coordinates  with  respect  to  the  curvilinear  coordinates.  With  the  for¬ 
mulae  of  Gauss  as  the  basis  for  the  proposed  elliptic  equations  it  can 
be  stated  with  certainty  that  every  smooth  surface  (second  order  dif¬ 
ferentiable)  must  satisfy  these  equations  for  any  allowable  coordinate 
system  introduced  in  the  surface.  Work  on  these  lines  was  initiated  by 
Warsi  [6-9],  and  an  exclusive  article  on  surface  grids  [10]  may  also  be 
consulted.  Various  other  contributions  by  using  the  proposed  equations 
have  been  made  [11-15].  An  important  point  to  be  made  here  is  that  all 
other  successful  elliptic  models  for  surface  grid  generation  are  a 
consequence  of  the  proposed  equations  (Eq.  (5)),  e.g.,  Garon  and 
Camarero  [16]  and  Thomas  [17]. 


The  proposed  equations  (Eqs.  (5))  can  be  used  to  generate  the 
Cartesian  coordinates  as  functions  of  the  curvilinear  coordinates  when 
the  surface  has  been  specified  either  analytically  or  by  discrete  data 
points.  Further,  the  boundary  lines  in  the  piece  of  a  surface  in  which 
the  coordinates  are  to  be  introduced  have  to  be  decided  in  advance. 
There  is,  however,  no  restriction  whether  the  piece  of  the  surface 
forms  a  simply-  or  multiply-connected  domain. 


For  arbitrary  shaped  bodies  in  which  only  discrete  data  points  of 
the  surface  are  given,  it  is  important  to  fit  a  global  equation  of  the 
form  F(x,y,z)-0  so  as  to  express  the  mean  curvature  of  the  surface  as  a 
function  of  x,y,z.  A  method  which  has  worked  for  many  cases  including 
the  case  of  a  fuselage  is  that  of  overlapping  piecewise  least-squares 
method.  Results  for  some  geometrically  difficult  cases  have  been  pre¬ 
sented  in  this  paper.  Also,  as  an  application  of  the  surface  grid  gen¬ 
eration  scheme  reported  here,  the  case  of  monoclinic  coordinates  has 
been  considered. 

2.  NOMENCLATURE 
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form  f6?B 


;  the  coefficients  of  the  second  fundamental 
the  surface  xv-const. 
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-  differential  operator  defined  in  Eq.  (3). 
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Jacobian  determinant. 

principal  curvatures  at  a  point  in  xv=const. 
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unit  normal  vector  on  xv=const. 
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rectangular  Cartesian  coordinates  x,y,z. 
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3D  curvilinear  coordinates;  i  -  1,2,3. 
2D  curvilinear  coordinates. 


coordinates  in  successive  transformations; 

0,1,2,... 
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Other  Rules: 

1 .  Repeated  lower  and  upper  indices  will  always  imply  summation  over 
the  range  of  index  values. 

2.  Greek  indices  except  v  take  values  (1,2)  or  (2,3)  or  (3,1),  while 
the  Latin  indices  take  values  1,2,3. 

3.  If  an  index  is  enclosed  in  parentheses  then  the  summation  conven¬ 
tion  (rule  1)  is  inapplicable. 
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3.  THE  MATHEMATICAL  MODEL 


The  mathematical  basis  of  the  present  formulation  along  with  the 
derivation  of  the  elliptic  equations  for  the  generation  of  surface 
grids  has  been  discussed  in  various  publications,  [6]  -  [10].  The  main 
steps  have,  however,  been  summarized  below. 

Usually  the  main  aim  of  grid  generation  is  to  generate  the  rectan¬ 
gular  Cartesian  coordinates  on  the  curvilinear  parametric  coordinates 
and  hence  equations  must  be  obtained  in  which  the  curvilinear  coor¬ 
dinates  appear  as  the  independent  variables.  Fortunately,  the  formulae 
of  Gauss  provide  a  direct  access  to  those  quantities  in  which  the  sur¬ 
face  curvilinear  coordinates  already  appear  as  the  independent  vari¬ 
ables.  If  xa  (a  assuming  only  two  values)  are  the  surface  coordinates 
in  the  surface  xv=const.,  then  the  formulae  of  Gauss  (cf.  [7,18])  are 


r  -  t6  r  ♦  n(v)b  , 
-,a6  a8L,6  -  ct  8 


(1) 


where  all  quantities  appearing  in  (1)  have  been  defined  in  Sect.  2. 
Inner  multiplication  of  Eq.  (1)  with  G  ga  results  in  the  vector  equa¬ 
tion  v 


Dr  +  G  (A^v)x5)r  . 
-  v  2  -,<$ 


(2) 


where 
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A^x6 
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For  the  purpose  of  imposing  a  control  on  the  distribution  of  coor¬ 
dinates,  we  set 
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(v) 
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aBp5 
S  oB' 


(4) 
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where  P^ft  are  six  arbitrarily  specified  control  functions.  For  a  thor¬ 
ough  difcussion  on  the  properties  of  P*  ,  and  also  on  the  r|lg- 
tionship  between  the  Beltramians  *  and  the  3D  Laplacians  V  x  , 
refer  to  Warsi  [10],  [19],  and  to  Sect.  5  of  this. paper. 

Substitution  of  Eq.  (4)  in  Eq.  (2)  yields  a  deterministic  set  of 
the  grid  generation  equations.  As  an  example,  if  x  =c=const.  is  the 
surface  in  which  one  wishes  to  introduce  the  coordinates  £  and  n,  then 
the  three  scalar  equations  from  Eq.  (1)  are 

Lx  =  X(3)R,  Ly  =  Y(3)R,  Lz  =  Z(3)R,  (5) 


where 


L  "  g223£5  ~  2gl23cn  +  gll3nn  +  P3S  +  QV 


P  "  g22P11  2g1 2P1 2  +  g11P22’ 


Q  =  g22P1  1  ~  2gl2P12  +  g1 1 P22" 


Here 


g11  =  g22/G3*  g12  =  ~gi2/G3'  g22  “  g11/G3* 


G3  =  g1 1 g22  (g12^  * 


(6) 


4.  NUMERICAL  SOLUTION 


Numerical  solution  of  Eqs.  (5)  has  been  obtained  by  using  point  and 
LSOR  iterative  schemes.  The  main  difference  between  the  grid  genera¬ 
tion  in  a  flat  space  and  in  a  curved  surface  is  due  to  the  appearance 
of  the  forcing  terms  on  the  right  hand  sides  of  Eqs.  (5).  The  quantity 
R  has  been  defined  in  Eq.  (3)  and  is  composed  of  the  principal  curva¬ 
tures  as  local  functions  of  coordinates  in  a  surface.  Thus  the  distin¬ 
guishing  features  of  the  body  geometry  are  contained  in  both  R  and  the 
unit  normal  vector  n.  Before  attempting  to  solve  Eqs.  (5)  it  is  impor¬ 
tant  to  feed  the  body-geometry  information  either  through  R  or  n.  Most 
of  the  computer  runs  for  different  shapes  reported  in  Refs.  [To]  and 
[20]  and  in  this  paper  have  been  conducted  by  prescribing  R  as  a  func¬ 
tion  of  x,y,z.  That  is,  once  the  equation  of  a  surface  either  in  the 
form  F(x,y,z)=0  or  z-f(x,y)  has  been  established,  the  principal  curva¬ 
ture  is  expressed  as  a  function  of  x,y,z  by  using  the  formula. 


kT  +  kTT  =  [(F2  +  F*)(2F  F  F  -  F2F  F2F  ) 
I  II  y  2  x  z  xz  z  xx  x  zz 


+  2F  F  (F  F  +  F  F  F  -  F  F  F  -  F  F  F  ) 
x  y  z  xy  x  y  zz  y  z  xz  x  z  yz 


+  (F2  +  F2)(2F  F  F  -  2F2F  -  F2F  )]/P3F2, 
x  z  y  z  yz  z  yy  y  zz  J  z 


(7) 


where 


2  2  2 

F  *  F +  F  . 
x  y  z 


As  has  been  noted  in  [10],  if  at  any  point  F  »  0,  then  the  other  form 
of  Eq.  (7)  is  obtained  by  interchanging  x  by  y,  y  by  z,  and  z  by  x. 
The  information  supplied  by  (7)  is  enough,  since  both  G~  and  n  are  cal¬ 
culated  by  the  iterative  scheme  itself. 

Much  of  the  present  authors’  efforts  have  been  directed  toward  ob¬ 
taining  a  numerical  scheme  which  can  be  used  for  arbitrarily  shaped 
surfaces  so  as  to  fit  a  function  F(x,y,z)=0  based  on  the  given  discrete 
surface  data.  Some  success  has  been  achieved  by  using  overlap- 
piecewise  least  squares  method.  In  each  piece  a  second  degree  polyno¬ 
mial  in  x,y,z  has  been  used.  Many  known  surfaces  of  second  degree, 
e.g. ,  ellipsoid,  hyperboloids  etc.  have  been  duplicated  by  using  this 
method.  This  method  has  also  been  used  to  duplicate  an  airplane  fusel¬ 
age  for  which  discrete  x,y,z  values  were  available.  Most  of  the  dis¬ 
tortion  appears  in  the  canopy  part  of  the  fuselage,  though  the  coor¬ 
dinates  generated  by  Eqs.  (5)  seem  to  be  smooth.  Figure  5  demonstrates 
the  slight  distortion  in  the  body  geometry  and  the  generated  grids  on 
the  fuselage.  Figures  1-5  show  the  capability  of  the  grid  generation 
equations  (Eq.  (5))  for  some  difficult  cases  in  each  of  which  the  ana¬ 
lytical  form  of  the  equation  F(x,y,z)=0  is  known  beforehand.  For  grids 
generated  on  other  known  surfaces  refer  to  [10]  and  [20]. 


5.  BELTRAMIAN  AND  LAPLACIAN  APPROACHES 


For  the  purpose  of  comparison  of  the  surface  grid  generation  model 
(Eq.  (5))  with  the  surface  grid  generation  by  using  the  Laplace/Poisson 
equations  in  3D  at  a  surface,  Warsi  in  [10]  and  [19]  has  established 
the  relevant  relationships  between  the  two  approaches.  This  analysis 
shows  that  if  the  transverse  coordinate  (the  coordinate  going  out  of 
the  surface)  is  orthogonal  to  the  surface  then  the  inverted  form  of  the 
3D  Laplace/Poisson  equations  and  Eqs.  (5)  are  the  same.  This  result 
does  not  demand  that  the  transverse  coordinate,  besides  being  orthogon¬ 
al  to  the  surface,  should  also  be  a  line  of  zero  curvature,  [17]. 

Let  5  and  n  be  any  parametric  coordinate  in  a  surface  with  5  being 
orthogonal  to  the  surface.  Then  as  shown  in  [10]  and  [19],  the  rela¬ 
tions  between  the  Beltramians  and  the  Laplacians  are  as  follows: 


a25 


v2e 


33 

?33' 


A2n 


V2n 


33, 

*33 


Further 


V2c  = 


.111  _ 

g33 


k(3)+k(3) 

I  II 


'S 


33 


and 


k(3)  +  k(3) 

I  II 


3G_ 


2/i^  C3  9C 


G3  g1  1g22  (g12) 


(8) 
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6.  COORDINATE  CONTROL 


The  specification  of  the  Beltramians  in  the  form  of  Eq.  (4)  implies 
that  the  control  functions  P^g  have  been  chosen  a' priori  by  the  user. 
In  all  of  the  computer  runs  shown  in  Figures  1-5,  we  have  taken  P^  =0 
for  all  the  relevant  index  values.  One  case  in  which  u 


=  0, 


2 

and  P22  is  a  specified  function  of  n  (refer  to  [10],  Eq.  (4.6))  is 
shown  in  Fig.  6  which  amply  demonstrates  the  role  of  the  control  func¬ 
tions. 


Looking  at  the  control  functions  as  coordinate  redistribution  func¬ 
tions,  it  has  been  shown  by  Warsi  [10],  [19]  that  for 
coordinate  systems  denoted  as  x?  _1 . 

expressed  in  terms  of  P°„,  , .  as  ‘ 

aB(m-l) 


two  successive 

and  X?  then  the  P°„.  ,  can  be 
(m)’  aB(m) 


p^  =  -t^  +rp°  +t°  l 

aB(m)  a6(m)  ep(m-l)  ep(m-1)J 


3x 


3x 


3X1; 


(m)  OA(m-1)  OA(m-1) 


eX(m-1)  °X(m)  dX(m) 


(9) 


where 


<5 

6(0) 


for  all  relevant  values  of  5,ct,8. 


7.  GENERATION  OF  NON-RECTANGULAR  COORDINATES 


In  some  problems  it  is  desirable  to  generate  non-rectangular  cooi — 
dinates  on  curvilinear  coordinates  in  the  same  surface.  The  non- 
rectangular  coordinates  can  be  spherical,  cylinderical  or  any  general 
coordinates  in  the  surface.  The  formulation  of  the  pertinent  equations 
follows  directly  in  a  simple  way  from  Eqs.  (5)  and  has  been  fully  de¬ 
scribed  in  [1],  p.  248.  Refer  also  to  [13]-[14]. 

The  equations  in  this  case  can  concisely  be  written  as 

Lu  =  J^u  (10) 

where  u  =  (u,v)  and 

L  =  a3r,  -  2b3  +  c3  +  (P3r  +  Q3  ). 

5S  5n  nn  3  5  n 

For  an  expanded  form  of  the  coefficient  a,b, ...  refer  to  [1].  Here 
(u,v)  are  the  non-rectangular  coordinates  while  £,n  are  the  desired 
coordinates.  Figure  7  demonstrates  the  use  of  Eq.  (10)  in  obtaining 
the  coordinates  on  the  surface  of  an  ellipsoid  forming  a  doubly- 
connected  region. 
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8.  MONOCLINIC  SURFACE-ORIENTED  COORDINATES 


A  monoclinic  coordinate  system  is  that  in  which  the  transverse  or 
the  outgoing  coordinate  from  the  surface  is  rectilinear  and  orthogonal 
to  the  surface;  the  other  two  coordinates  which  are  in  the  surface  can 
be  either  orthogonal  or  non-orthogonal.  The  capability  of  surface  grid 
generation  technique  discussed  in  Sections  3  and  4  coupled  with  the  oi — 
thogonal  and  rectilinear  nature  of  the  third  coordinate  provides  a  sim¬ 
ple  method  of  obtaining  3D  grids.  An  immediate  use  of  such  coordinates 
is  in  the  numerical  solution  of  the  Navier-Stokes  and  the  boundary  lay¬ 
er  equations.  Coordinate  systems  of  the  type  under  discussion  have 
been  considered  by  Hirchel  [21]  and  more  recently  by  Lee  [22].  The  es¬ 
sential  mathematics  is  simple  if  use  is  made  of  tensor  algebra. 

Let  a  coordinate  system  xct(x1=5,  x2=n)  has  already  been  generated 
in  a  given  surface  by  the  method  described  in  Sections  3  and  4.  We  now 
consider  a  neighboring  surface  parallel  to  the  original  surface  and 
denote  the  position  vector  of  any  point  of  the  parallel  surface  by  r. 
Then 


r  =  r(xa,5)  =  r(x“)  +  4»(?)n(xa). 


(ID 


where  5  is  the  outgoing  coordinate  from  the  original  surface  on  which 
5  =  0,  and  <(>  is  an  arbitrary  (user  specified)  function  of  5.  Using  the 
obvious  equations 


r  •  n  =  0,  n  •  n  =0,  n  •  n  =0, 

-,a  -  -  -,a  -  -,a 

we  conclude  that  n=n,  which  is  the  condition  of  parallel  surfaces. 
After  some  tensor  algebra,  the  geometric  quantities  of  the  parallel 
surface  are,  (cf.  [22]), 


*..8  ’  <’  '  K*  )g«B  '  12  -  <kI  +  kII)  +  !^8 


&„8  '  "  -  <kI  *  *  K*ecS' 


K  •  Ki 1  -  (kj  *  k  H  *  K*2)"', 


(12) 


where  an  overhead  bar  is  used  for  quantities  in  the  parallel  surface 
and  K  is  the  Gaussian  curvature. 

Using  the  preceding  algorithm  a  series  of  parallel  surfaces  can  be 
constructed  starting  from  the  data  of  the  given  surface.  Each  surface, 
starting  from  the  first  parallel  surface,  is  obtained  algebraically  by 


using  Eqs.  (12).  Lee  [22]  has  taken 


*u>  -  <£1>  kJ'JM 


where  J=1  corresponds  to  5=0,  JM=raaximum  5-steps,  and  k>1  is  a  chosen 
constant.  Figure  8  demonstrates  the  surfaces  for  a  piece  of  an  ellip¬ 
soidal  surface. 
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9.  CONCLUSIONS 


The  proposed  set  of  elliptic  equations  for  the  generation  of  sur¬ 
face  meshes  (Eqs.  (5))  have  been  solved  in  many  cases.  In  each  case 
the  resulting  coordinate  lines  are  sufficiently  smooth.  An  important 
aspect  of  the  algorithm  is  to  specify  the  equation  of  the  surface 
F(x,y,z)=0  based  on  the  discrete  input  data  for  a  surface.  This  func¬ 
tion  must  be  differentiable  at  least  to  the  second  order  for  it  to  be 
used  on  the  right  hand  side  terms  of  Eq.  (5).  Though  we  have  used  a 
least  squares  aproach,  our  research  in  this  area  is  still  continuing. 
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Figure  4, 


Coordinates  on  a  surface  having  a  saddle  point 
of  the  higher  type: 


xy(x2  -  y2),  -l<x<i,  -i<y<l 


Mi 


Figure  5.  Three  views  of  coordinates  on  a  fighter  plane 
fuselage. 


Figure  6.  Example  of  coordinate  contraction  in  a 
doubly-connected  region  on  the  surface  of  an 
ellipsoid. 
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Figure  7.  Coordinates  in  the  doubly-connected  region  on 
an  ellipsoid  by  using  Eqs.  (10). 
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Figure  8.  Monoclinic  coordinates.  J  *  1  is  the  basic 
ellipsoidal  surface  in  which  the  coordinates 
have  been  generated  through  using  Eqs.  (5), 
while  J  =  5,  9,  13  are  the  parallel  surfaces. 


I® 

\m 


Publications  under  Grant  AF0SR-85-01 43 


1.  "Numerical  Grid  Generation  in  Arbitrary  Surfaces  Through  a 
Second-Order  Differential-Geometric  Model",  by  Z.U.A.  Warsi, 
Accepted  in  J.  Comp.  Physics. 

2.  "A  Synopsis  of  Elliptic  PDE  Models  for  Grid  Generation",  by  Z.U.A. 
Warsi,  Accepted  in  Applied  Math,  and  Comp. 

3.  "Generation  of  Surface  Coordinates  by  Elliptic  Partial 

Differential  Equations",  by  Z.U.A.  Warsi,  NASA  Conf.  Publ.  2379, 
June  (1985). 

H.  "Numerical  Grid  Generation  Through  Second  Order 
Differential-Geometric  Models",  by  Z.U.A.  Warsi,  and  W.N.  Tiarn, 
11th  IMACS  World  Congress,  Oslo,  Aug.  (1985). 

5.  "Generation  of  Three-Dimensional  Coordinates  by  Using  the  Method 
of  Parallel  Surface",  by  J.H.  Lee,  MS  Thesis,  Mississippi  State 
University,  Dec.  (1985). 


18 


